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Abstract 

We construct asymptotic solutions of the functional Schrodinger 
equation for a scalar field in the Gaussian approximation at large 
proper time. These solutions describe the late proper time stages of 
the expansion of a meson gas with boost invariant boundary condi- 
tions. The relevance of these solutions for the formation of a disori- 
ented chiral condensate in ultra relativistic collisions is discussed. 
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Introduction. — Several recent studies have been devoted to the cal- 
culation of the evolution of a self interacting scalar field. As an example 
the evolution of a scalar field is an important ingredient in most inflationary 
models [p], ||, [3], £|, |[. Also classical as well as quantum calculations of the 
evolution of a scalar field have been considered in the context of the possi- 
ble formation of a disoriented chiral condensate in ultrarelativistic nuclear 
collisions §, 0, 1, 1, [10 . 



In this note we investigate the asymptotic forms of the quasi-one- dimensional 
quantum evolution of the meson fields which possess a symmetry with re- 
spect to the Lorentz boost along z-direction. Such a symmetry was orig- 
inally introduced by Bjorken|[TT|| in a hydrodynamical model of ultrarela- 



tivistic nucleus-nucleus collision based on the physical assumption of central 
plateau formation in the final particle distribution in the rapidity space. 

We take a simple model for mesons described by a self-interacting scalar 
field ip{x) and adopt the functional Schrodinger picture with the Hamiltonian 
density: 

w(x) = ^ 2 (x) + i(\Mx)) 2 + ^ 2 ( X ) + A^ (x) 

We use the framework of the Gaussian or mean field approximation, i.e. we 
assume that at each time t the wave functional describing the system is a 
gaussian type coherent state 

\& [<^(x)] = Mexp (i < ti\lp — <p >) 

x exp (- < (p - <p\{^q + iE)\(p ~<P >) > 

where G, S, <p, n define respectively the real and imaginary part of the kernel 
of the Gaussian and its average position and momentum. We impose Lorentz 
invariant boundary conditions i.e. we assume that observables depend only 
on the proper time variable r = y/t 2 — z 2 . We examine the behaviour of the 
meson wave functional for large values of the proper time variable. We will 
show that this limit can be worked out analytically by solving the linearized 
evolution equations. 

For more simplicity we restrict ourselves to the case of a single scalar 
field with an initial condition which reduces to a single pure state (not a 
statistical mixture). As will be shown below this case is already instructive 



2 



and general enough to provide a transparent interpretation of some of the 



numerical calculations performed recently by Cooper et al jl2| in the more 
general case of a sigma model with an initial state which is a statistical 
mixture. A more specific and more detailed comparison adapted to this case 
will be presented in a forthcoming publication. 



Mean Field Equations. — 

approximation read 0, U f| 

G = 2(G£ + £G), 



The evolution equations in the Gaussian 



2S 2 -i{-A + m^ + ^ 2 + ^(x,x)} 



-7T, 



tt = {-A + mg + ^> 2 + ^G(x,x)}£. 

o 2 



In these equations m is the bare mass of the scalar field and A the bare 
strength of the self coupling term which we write as A<y? 4 in the Lagrangian 
density. The vacuum state corresponds to the static solution 



0. 



7T 



0, <p = <p , G= l/2J-A + fx 2 , 



where the quantity \i has to satisfy the so- called gap equation [[Lj 

2 2 I 1 I ^2 

u = m n H — < x — , = x > +— o9 n . 

P 2 '2V-A + /Z 141 2^° 



(1) 



The previous equation requires a regularization scheme such as a discretiza- 
tion of the Laplacian operator on a lattice with a mesh size Ax = 1/A or a 
cutoff A in momentum space. To make the gap equation finite when the scale 
A goes to oo a popular prescription is to send the bare coupling constant to 



zero according to the formula [H 

1 



2A 



R 



1 1 

A + 16^ 



log 



2A, 
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This prescription is however not free from difficulties, see 0, [T3|, [14| . 

The mean field evolution equations can be written in a more compact 
form which has the further advantage of being manifestly covariant 

m 2 (x) = ml + ^ 2 (x) + 2" < X \S\ X > j 

{ n + mo + ^ 2 + ^ < x\S\x >}<p = 0, 
o 2 

where x = (xo,xi,X2,Xs) and where S is the Feynman propagator in the 
presence of an x dependent mass 

S 



□ + m 2 (x) + ie 



A proof of this equation can be constructed by using the functional methods 
of Cornwall, Jackiw and Tomboulis [[HJ . The kernel G of the Gaussian wave 
functional is related to the propagator S via the relation 

< x|G(t)|x >=< x\S\x > . 

To obtain the asymptotic form of the state of the meson field at large 
proper times we assume that the propagator S is close to its vacuum value. 
In this case we linerize the equations of motion around the vacuum value by 
writing 

S = — — 2(~\ i • = S o + S o i 5 ™ 2 ) S + ... (2) 
Introducing the Fourier decomposition 

5 m 2 (x) = [ — -A— exp(iq ■ x)5m 2 (q), 



(2vr) 4 

" = 1 (2^) 
equation (0) leads to 

< x\5S\x >= J — — exp(iq ■ x)Il (q 2 )5m 2 (q) 



f d A q 

5 <f(x) = (p(x) -<p = ——- exp(iq ■ x)5(p(q), 
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where ITo(g 2 ) is the standard polarization tensor 

n (g 2 ) = - J -^^So{p)S (p + q) 

1 n / 4A \ <f 
lbir z e A [i l fi 

The linearized equations of motion are thus found to be 

5 m 2 (q) = \(p 6(p(q) + ^Tl (q 2 )5m 2 (q), 
( -q 2 + /i 2 )^(g) + ^ U (q 2 )5m 2 (q) = 0. 



These equations are equivalent to those derived by Kerman and Lin 14 



in 



the functional Schrodinger picture. A compact and elegant construction is 
provided by our manifestly covariant formalism. 

When the momentum cutoff A goes to oo the quantity A n has a finite 
limit while the product A <po goes to zero and the first equation becomes 



^n (g 2 ) 



Sm 2 (q) = 0. 



The linearized evolution equations have two types of solutions. The first type 
of solutions are given by 

5m 2 = 0, 5 V (q) = f(q)5(q 2 -^ 2 ), 

where \i is the self consistent solution of the gap equation (|I|) and f(q) is 
a regular function of q = (qo,qi,q2,q3) to be specified later by boundary 
conditions. In the second type 

5m 2 (q) = (q 2 - ^ 2 )5<p dip = f(q)5(q 2 - M 2 ), 

One important difference with the solutions of the first type is that now M 2 
is a solution of the equation 

i - ^n (M 2 ) = o. 
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This equation is just the lowest order Bethe-Salpeter equation for the ip 4 field 
theory. For small enough values of the renormalized coupling it has a single 
solution M 2 i.e. a single bound state |T3J . 



To construct Lorentz boost invariant solutions inside the forward light 
cone (xq — x\ = t 2 — z 2 > 0, x = t > 0), it is convenient to introduce the 
light cone coordinates by 

Xq = r cosh 77, X3 = t sinhr], , 
q = acoshy, q 3 = asmhy, 



where a = ±W q" 2 — q 2 and y is the usual rapidity variable. Then the Fourier 
integral for dip for the first type of solutions is written as 

] /*oo /*oo P 

x exp [iar cosh(y — rf) — zqj_ ■ xj_] , 

where xj_ = (xi, x 2 ) and q_i_ = (qi, q 2 ). We seek the form of f(q) which gives 
(fi(x) independent of r\ and xj_. By noting that q 2 = a 2 — q^ does not depend 
on y, one can show by inspection that the required form is given by 

f(q) = 5 {2) (c l± )g(a) 

Indeed, the Fourier integral can now be performed explicitly and the result 
can be written by a linear combination of the Bessel function and the Neu- 
mann function of order zero, Jo an d N : 

(p(x) = y? + aJo(fiT) + PN (ht), 



where we have used the following integral expressions for J and iVo 16 



2 f°° 

Jo(x) = — rftsin(xcosht) 

TT JO 

2 r°° 

N Q (x) = / dt cos (x cosh t) 

ix Jo 

The two coefficients a and [3 are given by the value of the function g(a) at 
a = ±/i 

1 % 

« = \g{-f) - g(v)] , ^ = 32^3 [ffW + 



These two coefficients may be alternatively given by the initial conditions of 
<p{t) and its derivative at r = tq when the evolution of the meson field starts. 

The appearance of the Bessel (and Neumann) function may be easily 
understood if we rewrite the Klein-Gordon equation (□ + ji 2 ) ip(x) = using 



the light cone variables [12 



d 2 Id Id 2 d 2 2 
dr 2 r dr t 2 drf dx\ ^ 



ip(x) = 0. 



For (p(x) independent of rj and xj_, this equation becomes a form of Bessel's 
differential equations whose two independent (real) solutions are given by Jo 
and N . 

Boost invariant solutions of the second type can be obtained in a similar 
way: 

y? (x) = <p + aJ (Mr) + (3N (Mt), 

m 2 (x) = /i 2 + (M 2 - /j 2 ) [aJ (Mr) + (3N (Mt)} . 

We thus learn that for large values of the proper time the system returns, 
as postulated above, to the vacuum state for both types of solutions, which 
justifies a posteriori the approximation scheme we have developped. 

Discussion. — The first type of solution leads to a momentum distribu- 
tion of the form 1/y/k 2 + /i 2 which is time independent and boost invariant. 
It describes the Lorentz invariant expanding coherent state of a meson gas. 
The second type of solution involves the solution of the lowest order Bethe- 
Salpeter equation and can thus be interpreted as an expanding coherent state 
of two-meson bound states. Solutions exist only for a negative bare coupling 



constant, as discussed in reference 14 



Although our formulae were derived for the specific case of a single scalar 
field with pure Gaussian states and a particular renormalization scheme in- 
volving a negative bare coupling, they exhibit general features of boost invari- 
ant boundary conditions which are expected to hold for other more general 
models. In this respect it is instructive to compare our formulae with the 
results obtained numerically by Cooper et al ||12|| . These calculations con- 
cern the case of a sigma model with statistical mixtures as initial conditions. 
They are expected to be relevant to discuss the possible formation of a chiral 



7 



condensate in ultra relativistic collisions. They use a finite momentum cutoff 
and a positive bare coupling constant, a case to which our formulae are not 
a priori immediately applicable. There is nevertheless a striking agreement 
between our formulae involving Bessel functions of order zero and figure 5 
of Cooper et al which shows the variation of the square mass /x 2 + 5m 2 as 
a function of the proper time r. It is interesting to note that the results 
of Cooper et al are well described by a single Bessel function whereas our 
previous analysis suggests a superposition of two Bessel functions with char- 
acteristic scales /i and M. One should remember however that since Cooper 
et al have a positive coupling constant there are no two- meson bound states. 
This explains why the scale M is absent but still leaves one unanswered ques- 
tion. Indeed without bound states our formulae predict no variation at all 
in the square mass. The answer to this question is that our derivation was 
made in the case of an infinite momentum cutoff whereas Cooper et al use a 
finite one. This difference produces the non vanishing value of 5m 2 . 

By looking at the results of Cooper et al it is also worthwhile to note that 
some of the initial conditions make the small amplitude approximation valid 
for the whole evolution. Our formulae also explain nicely why increasing 
the formation time leads to larger asymptotic oscillations as a result of the 
1 / yfr behavior of the Bessel functions. They do not explain why there are 
in this case less instabilities since this involves the small proper time domain 
in which our approximation scheme is expected to break down. Our analytic 
asymptotic expressions appear nevertheless as a useful reference to analyze 
some of the physics of the boost invariant expansion of a meson gas. Results 
using the sigma model with a finite cutoff will be presented in a forthcoming 
publication. 
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